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Abstract 

Some general properties of perturbed (rational) CFT in the background metric of sym- 
metric 2D sphere of radius R are discussed, including conformal perturbation theory for the 
partition function and the large R asymptotic. The truncated conformal space scheme is 
adopted to treat numerically perturbed rational CFT's in the spherical background. Numer- 
ical results obtained for the scaling Lee- Yang model lead to the conclusion that the partition 
function is an entire function of the coupling constant. Exploiting this analytic structure 
we are able to describe rather precisely the "experimental" truncated space data, including 
even the large R behavior, starting only with the CFT information and few first terms of 
conformal perturbation theory. 

1. Motivations 



There are several reasons to be interested in euclidean 2D field theory in spherical back- 
ground. Apart from conceptual interest in relativistic field theory on curved spaces, the 
particular example of sphere presents an apparent simplification of a symmetric space where 
physics is independent on the space position while being different from that on the infinite 
flat plane. 

It has been proven useful to study properties of interacting field theory through finite-size 
effects (see e.g. [jT], 0, [J). The compact geometry of sphere naturally provides a settlement of 
finite size problem, different from the Casimir effect which has been mainly studied so far. 
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Exact constructions of various conformal field theories (CFT) in 2D gave rise to a very 
fruitful idea to consider non-conformal (both massive and massless) field theory models as 
certain perturbations of CFT's by relevant operators 0j. One considers formal actions like 

A x = Acft + 77^/ <t>{x)d 2 x (1.1) 
lix J 

where <fr is some relevant primary field in the theory of dimension A and A is the corresponding 
coupling constant. Even if such perturbation leads to a field theory which is not integrable, 
the perturbation theory in A, called the conformal perturbation theory (CPT), may give 
valuable information about the short distance behavior of the model. One of the main 
problems of CPT is related to infrared divergences of the perturbative integrals, which result 
in certain non-analyticity in A. Finite geometry of a sphere imposes an infrared cutoff. It is 
natural to expect that on a finite sphere observables are analytic in A at A = 0. Below I'll try 
to argue that in certain cases the CPT series on a finite sphere is absolutely convergent and 
therefore results in entire function of A. Although there are no obvious reasons to expect any 
sort of integrability on the sphere even if the model ( |1.1| ) is integrable on the flat plane, the 
analytic property suggested is extremely restrictive. As we'll see before long, it permits to 
get a lot of information starting from few first CPT coefficients and some general information 
about A — > oo asymptotic. 

There is an approximate way to treat the actions like (|1 . 1|) called the truncated conformal 
space (TCS) approach ||. Sometimes it gives rather precise numerical information even 
about the large-distance properties of the field theory ||. Spherical geometry seems to be a 
very natural environment for the TCS approach. 

The main motivation, at least for me, to study perturbed CFT's on fixed sphere is the 
long standing challenge of 2D quantum gravity. There are all reasons to believe that the field 
theoretic approach based on the Liouville field theory is a relevant way to understand 2D 
gravity at least in the so-called "week coupling regime" . However, up to now only very basic 
information about general scale dimensions and simplest correlation functions is reached by 
the field theoretic means, even in the simplest case of spherical topology. This has to be 
compared with very reach "experimental" results coming from the matrix models and similar 
"discrete geometry" approaches (see e.g. || for most recent achievements in the field). I hope 
that better understanding of the field theory on a fixed symmetric sphere can be a starting 
point to approach the problem at least in the semiclassical limit where the central charge of 
the massless matter modes is large negative and the geometry becomes "rigid" . 

2. CFT in curved background 

The most advanced understanding of the 2D field theory in curved background g a b{x) is 
reached in conformal field theory. This is due to the following two simple properties postu- 
lated as a basis of CFT construction 

Stress tensor anomaly. In 2D CFT the trace of the stress tensor T a b(x) is in fact a 
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c-number (i.e., proportional to the identity operator) and reads explicitly 



6{x)=g ab T ab = -^R (2.1) 

where R is the scalar curvature of the background metric and c is certain (real) number, 
called the central charge. The central charge is the basic characteristic for any particular 
CFT. 

Primary fields. Among the set of local fields {$} in the theory there is a subset of so 
called primary fields <pi which transform very simply (here the primary fields are supposed 
to be scalars) 

54>i{x) = -^<j>i(x)5ip{x) (2.2) 
under the Weil variations of the metric 

Sg a b(x) = g a b(x)S<f(x) (2.3) 



In fl2.2p Aj are (real) numbers characteristic for the primary fields, called the scale dimensions 
(or conformal weights). All other local fields can be found in the operator product expansions 
of primaries and the nontrivial stress tensor components. Especially simple are the so called 
rational CFT's which involve only finite number of primary fields (and therefore contain 
finite spectrum of dimensions Aj). 

Conceptually these two properties (which were in fact abstracted from explicit calcula- 
tions in certain simple examples like free field theories) are enough to develop the whole 
structure as rich as the conformal field theory. 

In 2D one can always choose locally an "isothermal" coordinate system where the metric 
is "conformally flat" 

9ab(x) = e^6 ab (2.4) 
In the case of sphere such coordinates (z, z) can be chosen almost globally, covering all the 



surface except for the "south pole" z = oo. Simple transformation low (2^2) allows to express 
any correlation function of primary fields through the correlation functions on infinite flat 
plane with g ab = 5 ab (again all fields are implied scalar) 

n 

(Mxi) ■ ■ ■ = J] e ' AMXl) ■ ■ ■ <PnM) v=0 (2.5) 

i=l 

In this paper we'll concentrate on the metric of sphere of radius R. The Weil factor 
exp (p(z, z) in ( |2.4j) can be chosen in the form 

= (2.6) 

1 + zz Y 
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The scalar curvature is a positive constant 

R = -Ae^ddif = 2R~ 2 (2.7) 

Metric ( |2.6|) is invariant under the 0(3) group of coordinate transformations 

az + b 



a — bz 



(2.i 



where it is convenient to set aa + bb = 1. 

The stress tensor anomaly ( J2.1| ) readily prescribes the R dependence of the partition 
function 

Z CFT (R) = R C/3 Z (2.9) 

Here Z = ^cft(I) is the CFT partition function at R — 1. As long as I know for this 
multiplier any theoretical predictions are lacking. In what follows I'll consider only the 
ratio Zqft(R) / Zo. A simple consequence of eq. fl2.3|) is the R dependence of any correlation 
function of primary fields on the sphere 

(M*l) ■ ■ ■ M*n)) R = R~ 2E7 ^ Al (faM ■ ■ ■ MXn)) R=1 (2.10) 

3. Perturbation 

In non-trivial background (O) perturbed action fll.ip is replaced by 



Ax = Acft + ^- [ <P(x)e^d 2 x (3.1) 
2tt J 



where <fi is again a relevant primary field of dimension A < 1. 

Formal expression (|3.1| ) gives rise to perturbative expansion in A. Consider e.g., the 
perturbed partition function Z\[ip\. Formally 

ZM . ^ (-A) 



00 (—\\ k r n 

s / • • • ^ xn)) - n ( 3 - 2 ) 

i— n v / ' " i— i 



where the correlation functions are computed in unperturbed CFT but in non-trivial back- 
ground ip. With the use of (|2.5|) one can "flatten" these correlation functions, reducing ( |3.2j ) 
to an expression which involves only the CFT correlation in trivial flat background ip — 



00 (—\) k f n 

£ / ■ ■ ■ *M>gft II e^^d^ (3.3) 

k=0 \ ' ' i=l 



ZcFTi^p 

On the sphere ( |2.6p the i? dependence is summarized as follows 

Zx(R) 



ZoRc/3 - 4h) (3.4) 
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where the dimensionless combination 

h = \(2R) 2 ' 2A (3.5) 
is introduced. Eg. ( |3.3|) gives z(h) as a regular power series in h 



n=0 



where explicitly 
z = 1 



1 /" n 

(0(^0 . . . 0«i) CFT n (1 + xx v~ 2A (3 ' 7) 



(2vr) n n! 



7T / T-r (11/ 



\2-2A 



at n > 



In the last equation the symmetry (|2.8| ) is used to eliminate one of the integrations. In CFT 
two- and three-point functions are fixed up to normalization constants. Assuming the usual 
in CFT normalization of primary fields we have 

(0O))cFT = 

(0Ol)0(>2)>cPT = 1 [4A ( 3 - 8 ) 

Fi - x 2 \ 



where is the corresponding structure constant, which is explicitly known in exact CFT 
constructions. First two coefficients in (|3.6| ) are readily calculated 



zx = (3.9) 
1 



Z-2 



8(1 -2A) 

while for z 3 we have 

= r d 2 y 1 d 2 y 2 

Zz 6(2nfJ [(1 + ym)(l + y 2 y 2 )} 2 - 2A [yiyiy2y2(yi - V2)(yi - V2)] A l ' 1 



C w r(2-3A)r(2- A) 
24 r(4 - 4A) 



v' A dvF(A, A, 1, u)F(2 - 2A, 2 - 3A, 4 - 4A, 1 - v) 



Higher perturbative integrals in ( |3.7|) involve four and more point CFT correlation functions. 
These are not that universal as (|3.8j ) and require more detailed information about a particular 
CFT model. 
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If the series ( |3.6| ) is convergent, it defines the perturbed partition function as a function of 
R. In what follows we'll see an evidence that in certain cases this series is indeed convergent 
and, moreover, absolutely convergent, so that z(h) is an entire function of h. 

In any case the perturbative development (|3.6| ) describes the R — > asymptotic of the 
partition function Z\(R)/Z . Let's now turn to the opposite R — > oo, or h — ► oo limit. 



4. Large R asymptotic 

It seems intuitively clear that when R ^> m _1 (where m ~ A 1 '^ 2-2 ^ is the typical mass scale 
in the perturbed model) the local physics is almost the same as in infinite flat space-time. 
In particular the leading exponential asymptotic of the partition function Z\(R) is expected 
to be governed by the specific vacuum energy £ vac in flat space 

log Z X (R) ~-4nR 2 £ yac + ... (4.1) 

This vacuum energy is an important non-perturbative characteristic of a field theory model. 
For dimensional reasons 

£ vac = -A\ 1/(1 ~ A) (4.2) 

where A is a dimensionless number. In integrable field theory this number is typically known 
exactly (see e.g., the example below). 

To get an idea about further R —>■ oo corrections let's start with the relation 

and take the rescaled metric 

As R — > oo this Weil factor is trivial e v = 1 and (9) coincides with the expectation value on 
infinite flat plane 

Wflat = 4 ^vac (4.5) 

reproducing the leading asymptotic ( |4.1| ). For next corrections we should take into account 
the stress tensor anomaly (2.1) 



9 = ~ + fla t (4.6) 

as well as the deviation of the metric from flat around the location of the operator 9{x). 
Let's put x = so that 

w -i z z 3(zz) 2 , A _. 
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General variation formula (X is any observable) 

5(X) = -— [ (9(x)X)e^ x) 5^(x)d 2 x (4.8) 
4tt J 

to the first order in 1/R 2 gives rise to the following R~ 2 correction to (9) 

(^(0)) sphcrc = <0(O)) flat + ^55 / <^M0)> flat * V* - + 0(7T 4 ) (4.9) 



/sphere \ u V";/flat r / \ u V^; v V";/flat ^ " qj^ 2 

The famous c-theorem sum rule |7j allows to evaluate the integral 

A-7CC 

(e(x)e(o)) aat x 2 d 2 x = — (4.io) 

The anomalous R~ 2 term is cancelled out and the IR corrections to (9) appear at the order 
.R -4 , so that 

b 2b 

(^(0)) S phere ~ ^™ + ^ + j£ + ■ ■ ■ C*-") 

with some (dimensional) coefficients bi, b 2 , ....Thus, as it can be intuitively expected, at 
R — > 00 there is no power correction to Z\(R) and the IR expansion has the form 

log ~ -^R 2 £^ + \og( Zoo ) + !jL + + . . . (4.12) 

Note that the sum rule Q4.10|) plays essential role in the cancellation of the log R correction. 
The (dimensional) integration constant 

Zoo = \-^-^ Zoo (4.13) 

where is a dimensionless number. Presently I don't understand neither physical meaning 
of this constant nor how it can be predicted on the basis of flat field theory. Nevertheless 
this parameter seems to be an important characteristic of the perturbed CFT on a sphere. 
In terms of the variable ( |3.5| ) asymptotic ( 4.12|) reads 

logz(h) = 7rA/iV(i-A) + l g(2 c / 3 z 0O ) °— log h + ^h^ 1 ^ + a 2 h~ 2 '^ + ... 

6 — 6A 

(4.14) 

where A is defined in (p|) and a x = 4M 1/(1_A) , a 2 = 166 2 A 2/{1_A) etc. 

Unlike z^, coefficients a 1; a 2 , etc., (or bi,b 2 , ... in ( [4.12| )) can be expressed in terms of 
higher correlation functions of 9 in flat space-time. For example 

h = J {e{x l )e{x 2 )em^ t x 2 l x 2 2 d 2 x l d 2 x 2 (e(x)9(o)) 6&t x*d 2 x (4.15) 

In an integrable model the flat space correlation functions in principle can be constructed in 
terms of the form-factors of the operator 9, which are typically known exactly. However, even 
numerical calculation of the three-point function (or multipoint ones for the next coefficients) 
in this framework presents a separate problem. 

Anyway, in the lack of any exact non-perturbative approach, a kind of "experimental" 
tool to measure the observables under consideration is very desirable. 
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5. Schrodinger picture in sphere geometry 



In [[| the truncated conformal space (TCS) approach has been used to evaluate numerically 
certain characteristics of 2D perturbed CFT's. In this paper and a number of subsequent 
works || it was demonstrated that TCS is reasonably effective for many 2D models. In this 
section a way to apply similar approach in the spherical geometry is discussed. 

For our purpose it is convenient to map the projective coordinates (z, z) on the sphere 
to the "cylindric" ones (t, a) 

z = exp £ = exp(t + icr) (5.1) 



where the sphere metric (|2.6| ) has the form 



e <p(t,*) = (5 _2) 
cosh t 

Consider the space of states on a circle of fixed a. We will be interested in the "time" t 
evolution of these states and, in the TCS spirit, choose as the (time independent) basis the 
set of CFT states \a) (here index a runs over all the CFT states, primary and descendent). 
For simplicity we suppose that the basis is orthonormal. Any state (Schrodinger picture is 
implied where the states are time dependent) 

|¥(f)> = £V(f) |a) (5.3) 

is described thus by the wave function ty a (t). Forgetting for the moment the conformal 
anomaly (whose effect on the partition function is summarized to the multiplier (2i?) c / 3 in 
( p.9| )) we consider the following time dependent Hamiltonian 

D(t) = D Q + g(t)V (5.4) 

which generates translations in the t direction. Here 

D = L + L (5.5) 

and Lo, Lq are the standard Virasoro generators acting on the CFT states. The interaction 
part of the evolution is constructed as 



V = — <j)(0,a)da (5.6) 
2vr Jo 

The time-dependent coupling constant g(t) in ( PT4] ) reads in terms of the "effective coupling" 
TBI) 

9{t) = - (2coshV 2A (5 - T) 



S 



Matrix elements V \a) = J2 b B b a \b) between CFT states are essentially the CFT structure 
constants 



B b a = ^ (f\<f>(0,a)\i)da = Cl a (5.8) 
^ Jo 

They are constructed explicitly in solvable CFT models. The t-evolution of the wave function 
is described by the time-dependent Schrodinger equation 

~ Mt)) = D(t) |*(f)> (5.9) 
In the CFT basis we arrive at the infinite dimensional system of linear differential equations 

b 

Different solutions to this system correspond to various operators placed at the points 
t = — oo and t = oo (i.e., at the "north" and "south" poles of the sphere). In particular, the 
solution ^vac(^) determined by the initial condition 

<C c (-oo) = 5? (5.11) 

where a — I is the CFT state corresponding to the identity operator, describes the state 
radiated by smooth (no field) north pole. The reduced partition function ( |3.4| ) reads from 
this solution 



z{h) = lim *Ut) (5.12) 



t^oo 



It is easy to verify that the formal development (which implies that all infinite sums over 
intermediate states are convergent) of the solution \l/ vac in powers of h results in the same 
perturbative series ( |3.6| ), (|3.7|) . From other components of the wave function one can read-off 
the one-point functions at the south pole 

^M = ^*Jp£M ,5,3) 

Here $ a is a scalar field corresponding to CFT state \a) and A a is its dimension. Of course, 
this limit is in general divergent and requires standard UV renormalization, necessary to 
define the perturbed fields (see e.g. [|, §). This point will be discussed in more details in 
rOfl while here we'll concentrate only on the partition function, where the limit ( |5.12p is well 



defined. 

The TCS idea is quite simple: truncate the infinite dimensional CFT space of states 
up to certain maximal dimension and then treat the resulting finite-dimensional problem 
( |5.10| ) numerically. Previous TCS experience shows that such procedure is often convergent 
(sometimes rather fast) with the increase of the truncation dimension. 
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6. Scaling Lee- Yang model 



Scaling Lee- Yang model (SLYM) [|ll[j is often used as a testing tool for various approximate 
approaches in 2D field theory. It is probably the simplest (excluding, of course, the free field 
theories) example of perturbed CFT. The model arises as a perturbation of the non-unitary 
CFT minimal model Ai(2/5). This CFT is rational and contains only two primary fields, 
the identity / of dimension and the basic field <fi = $1,3 of dimension A = —1/5, the central 
charge being negative c = —22/5. The basic operator product expansion reads 

0(x)0(O) = (xx) 2/5 I + (xx) 1/5 C^0(O) + . . . (6.1) 

where the three-0 structure constant is purely imaginary = in and 

1/2 

r 2 (i/5) 

- 1.911312699. . . (6.2) 




5r(3/5)r(4/5) 



^slym = ^cft + — / (t>{x)e^ x) d 2 x (6.3) 

is what is called the scaling Lee- Yang modelQ. In flat space-time (<p = 0) it is integrable. 
Its particle content and factorized scattering theory are known exactly |~2|. The spectrum 
contains only one massive particle of mass m. 

Integrability allows to obtain many exact results about SLYM. E.g., the relation 



m = k\ 5/12 



between the mass m and the coupling constant A is known |13 

r(4/5)r(3/5) i 5/24 



k ^ 

" r(5/6)r(2/3) [ioor(i/5)r(2/5)_ 

The bulk vacuum energy is also found exactly P, IT 



1.2288903248 . . . (6.4) 



m 



£- = -J^ (6.5) 



so that the dimensionless parameter A in (|4.2|) is 



A exact = ^= = 0.2179745... (6.6) 



4 Notice that the coupling A in ( |6.3| ) is taken purely imaginary. This is nessessery to make the perturbed 
theory real. 



10 



7. TCS at low truncations 



In this section TCS for the Schrodinger equation (|5.10| ) is applied to the scaling Lee- Yang 
model. This model is especially convenient for TCS thanks to the small number of pri- 
mary fields and high degeneracy of modules at low levels. This results in low-dimensional 
truncated spaces even for relatively large truncation levels and makes the finite dimensional 
problem especially simple for numerical treatment. E.g., truncation level 5 (which means the 
maximum scale dimension in the truncated space 2A max = 10) gives 17-dimensional space 
of states. 

Up to level 5 the matrix elements B{ can be found in ||. Truncated system (|5.10 ) was 
solved numerically and the reduced partition function z(h) is evaluated as ( |5.12| ). In fig.l 
the combination 



Zx(R) 
Z 



exp(4 7 r J R 2 ^ vac ) = ^ c (oo)R- 22/15 ex V (4nR 2 £ v 



(7.1) 



is plotted as a function of dimensionless variable mR = (k/2)h 5 ^ 12 , where m is the mass of 
SLYM particle. Data for different truncation levels are presented to illustrate the conver- 
gence. 




-0.6 



-0.8 



Figure 1: Combination (|7.1|) for truncation levels 5 (full curve) and 4 (small circles). The 
IR fit with formula ( |7.6|) is plotted in crosses. 
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First interesting feature we observe is a zero of the partition function at real positive 
R = Rq. Its position is estimated as 

rn.Ro = 0.88964 ... or h = 2.43083 . . . (7.2) 

At small R the data are well fitted by the perturbative series fl3.6| ) 

Zx{R) ~ z(h) = l-^-h 2 - z 3 h 3 + ... (7.3) 



Z R~ 22 / 15 K ' 56 
where Z3 can be computed numerically 

Z3 = ^^r^/S) 7 ^ jf vl/ * dvF (-V5, -1/5, 1, v)F(12/5, 13/5, 24/5, 1 - v) (7.4) 
= 0.0229899 . . . 

In the relatively stable interval 2.0 < mR < 3.0 the level 5 data were fitted, according to 
Q4.12j ), by the formula 

Z X (R) exp(47ri? 2 £ vac ) f 1+ h\ (7 5) 



Z 00 v R 2 

The best fit is achieved at 

Zoo = (-0.92 ± 0.02)m" 22/15 (7.6) 
h = (-0.72±0.2)m" 2 

For dimensionless constants and 61 in ( |4.14| ) this gives 

z^ = -1.25 ±0.03 (7.7) 
at = -1.9 ±0.5 

These numbers enter the h — > 00 asymptotic expansion 

log z{h) = uAh^ + log(2- 22 / 15 z 00 ) + 11 log /i + ai h- 5/6 + a 2 h- 10/e + ... (7.8) 

18 

where A is given by (|6.6[) . 

Analytic continuation for negative values of h is obtained by straightforward treatment of 
the truncated linear problem ( 5.1 0|) with negative h. The negative h data show an oscillating 
pattern presented in fig. 2, where the values of z(h)exp (tiA cos(7r/6)(— /i) 5//6 ) (— h)~ u / 18 are 
plotted (see below for the asymptotic at h — > —00) against —h. 

Numerical estimates for first few negative zeros h n , n = 1, 2, . . . are gathered in the first 
column of table |I[ It seems natural to expect infinite number of zeros on the negative h axis. 
Their positions are quite important for as we'll argue below z(h) is entire function of h and 
therefore essentially determined by the location of its zeros. 



12 



_•) 5 1 1 1 1 1 1 1 1 1 

'0 10 20 30 40 50 60 70 80 

Figure 2: The negative h data for z(h) exp ( it A cos(tt/ 6) (—h) 5 ^ 6 ) (— /i) -11 ' 18 at truncation 
levels 5 (continuous) and 4 (circles) 



8. Analytic considerations 



For any truncated finite dimensional linear problem ( |5.10| ), (|5.11 ) and (|5.12 ) the resulting 



reduced partition function z tTunc (h) is apparently an entire function of h. Let's suppose that 
this property holds also for the exact function z(h). We admit also that all zeros of z(h) are 
real, the first one ho being positive and the rest h n , n — 1, 2, . . . negative (and accumulating 
at h = — oo along the real axis). This implies in particular that the asymptotic ( |7.8|) holds 
in the whole complex /i-plane as \h\ — > oo, excluding the negative real axis arg/z = ±7r. This 
analytic structure combined with the perturbative information collected in ( |3.4| ) turns out 
to be quite restrictive. 

Asymptotic (|7.8|) allows to estimate the leading behavior of the locations of zeros h n at 
n — > oo 



1 



6/5 

-h n =[jln-^)+0(n- 1 )) (s.i; 



In table p] this asymptotic estimate is compared with the approximate TCS data for first 
several zeros. The asymptotic does surprisingly good even for the first negative zero h 1 . 
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n 




leading asymptotic 


sum rules 





2.43083 




2.43070 


1 


-11.762 


-12.41 


-11.7731 


2 


-30.439 


-30.66 


-30.2346 


3 


-50.60 


-51.05 




4 


-72.38 


-72.94 




5 




-95.98 




6 




-120.0 




7 




-144.9 





Table 1: Zeros of the partition function z(h) estimated by TCS, leading asymptotic ( |8 . 1| ) 
and sum rules (|8.3|) 



With (gy) the series 

n=0 n 

is convergent and defines the logarithmic derivative through the positions of zeros h n . From 
]3|) we have log z(h) = —T^h 2 — z^h 3 + . . . and therefore h n satisfy the following sum rules 



oo 1 



K 

n=0 " 

(8 ' 3) 

71=0 " 



= 32 3 = 0.0689697... 



n=0 ' ™ 



A quick inspection of table [T] shows that the leading asymptotic ( |8.1| ) agrees the TCS data 



with at least 1% accuracy even for the forth zero h 3 . Let's take this asymptotic expression 
as the exact one and use the sum rules to recalculate the first three zeros. This results in 
the numbers presented in the forth column of table [I]. They are impressively close to the 
"experimental" positions measured by TCS. This numerical observation can be considered 
as a strong support of the suggested analytic structure of z{h). 

Once all zeros are located with enough precision, the partition function can be restored 
as the convergent Weierstrass product 

"<*>=n(i-£)«p(£) ("> 



(due to (|8.1|) and (|8.3j ) the exponential multiplier here is not necessary and added to improve 
the convergence). The resulting function z(h) based on the above approximation for the 
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-1 ' 1 1 1 1 1 1 1 1 

01 2345678 

Figure 3: Weierstrass product (|S.4|) for Z\(R) exp(47ri? 2 £ vac )/Z based on the zeros in table 
|l] (full curve) against the level 5 TCS data. 



zeros h n is compared with level 5 TCS data for both positive and negative h in figs. 3 and 
4 respectively. 

In principle one could improve the precision using next-to-leading correction to the 
asymptotic ( |8.1|) 

- ft »KK""9 + M^ + 0( "" 2) ) 6/5 (8 ' 5) 

with "experimental" value ( |7.7[ ) for the coefficient a±. Here I prefer not to go along this 
line. One reason is the low precision in the estimate (|7.7f) . What is more important, the 
construction of z(h) via (|8.4|) haven't used any numerical TCS data at all, the letter being 
only the basis for certain hypotheses about the analytic structure. Moreover, the only 
non-perturbative information (which in fact comes from the integrability of SLYM in flat 
space-time) entering our above calculations is the exact vacuum energy ( |6.6| ). In the next 
section we'll see that even this income can be omitted (for the price of some tolerable loss of 
precision) and z(h) is restored rather accurately on the only ground of few first perturbative 
coefficients. 
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Figure 4: The same for negative values of the effective coupling constant h. Values of 
z(h)exp (7rAcos(7r/6)(-/i) 5/6 ) (-/i)" 11 / 18 are plotted. 

9. Integrability off 

Let's artificially switch off the exact information ( |6.6| ) and take this number as an undeter- 
minate parameter in the asymptotic expansion ( |7.8|) , together with and ai 5 a 2 , . . . . Of 
course, we continue to keep in mind the analytic structure of z(h) described above and de- 
note the zeros h n , n = 0, 1, 2, . . . as before. In particular the large n asymptotic (|8.1|) for h n 
still holds. 

To handle the zeros it is convenient to define (at Re s > 1) their zeta-function 



At 12/5 > Res > 1 it can be computed from the partition function z(h) as 

aY(s) = sin^/6) r°° h _ 5s/6 dlogz(h-i0) dh 
7T J dh 

where the contour of integration is shifted to agree the branch chosen in ( |9.1| ). Vise versa 




(9.1) 




(9.3) 
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where the integration contour goes to the left of the pole at s = 6/5 (in fact at s = 6/5 
there is no pole because Ciy(-s) vanishes at this point due to the sum rules, so that the first 
pole in the right half-plane appears at s — 12/5) and to the right of the pole at s — 1. The 
following properties of Cly(s) are readily figured out: 

1- Cly(~ Qn/5) = for n — 1, 2, 3, This is to avoid wrong powers of h in the asymptotic 

expansion (|Tlf). 

2. Cly(s) has simple pole at s = 1 with res s= i Cly(s) = A/2. This provides the correct 
leading asymptotic in (|7T8|) . 

3. Cly(O) = 11/18 to ensure correct log/i contribution in the asymptotic, while the con- 
stant term requires 

4- C[ Y (0) = ^log^oo-^log2. 
6 9 

5. Cly(s) has simple poles at all integer negative s = —n, n = 1,2,... and res s= _ n Cly(s) = 
7r _1 nsin(57rn/6)a n where aj, a2, . . . are the subleading coefficients in (|7.8|) . 

6. For all positive integer k = 1, 2, . . . we have by definition Cly(6/c/5) = Y^=o(~ h n )~ k , 
i.e., the numbers determined by the conformal perturbation theory. In particular (see 

(ED) 

Cly(6/5) = 
Cly(12/5) = 5/28 (9.4) 
Cly(18/5) = -323 = -0.0689697. . . 

The simplest approximation here is to replace all h n for n = 2,3,... by the leading 



asymptotic (8J.) while leaving A, h and hi as unknowns. In other words we approximate 
Cly(s) as 

^ = (4^ + ( 2 ) C(S ' 17/9) (9 ' 5) 

where ((s, a) = X]^=o( n a )~ s ^ s ^ ne usua l R-iemann zeta-function. This expression identi- 
cally satisfies property 3 while the sum rules ( ft.3p can be considered as a system of equations 
for undeterminate parameters A, ho and hi. There are two real solutions. The first one 

A = 0.218156. . . 

ho = 2.43068 . . . (9.6) 
hi = -11.7540... 

is really impressive when compared with the exact value ( |6.6|) for A and the "experimental" 
results for h n from table [I]. It's not clear if this remarkable precision is accidental. I'd like to 
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remind that it's not a first instance of miraculous numerical success while applying different 
approximations to this particular perturbed CFT || H . 

With these data in hand we can estimate also the constant evaluating Cly'(O) 

222/15TY17/OW5 

n - -1.22491... (9.7) 



/ i0 /ii(A/2) 5 / 3 (2vr) 3 / 5 
again in reasonable agreement with the TCS analysis. 



Approximation (|9.5| ) seems to be quite accurate in the right half-plane of s. However it 
has wrong analytic properties in the left half-plane. For example it has no pole at s = — 1 
(as well as at all s negative integer) and doesn't vanish at s = — 6n/5, n — 1,2, ... . E.g., 

/ 9 \ 6 / 5 

CS(-6/5) = -ho -h 1 + f - J C("6/5, 17/9) = -3.28339 . . . (9.8) 

One can try to cure the situation near s = — 1 introducing a pole at s = — 1 by hand, i.e., 
taking say (this ansatz is obviously inspired by eq.(|8.5|)) 

p 5iirs/6 i / A\ S n / A\ S+1 

C$00 = -pr + + (2 ) ^ 11 ^ - ^ s (2 ) C(s + 2 ' 17/9) (9 - 9) 

with a new undeterminate parameter a±. Then, with the previous values of A, ho and h\ we 
find that Cly(-6/5) = if 

01 w -2.060 (9.10) 

in tolerable agreement with ( |7.7| ). It can be verified that this correction doesn't violate 
essentially the sum rules ( |9.4|) while is slightly changed z^ = —1.26190 .... 

In principle one can continue the procedure and cure the value of Cly(— 12/5) by adding 
more terms and trying to correct the analytic structure near s = —2, s = —3 etc. Due to 
the asymptotic character of the expansion ( |7.8j ), such iterative procedure is hardly conver- 
gent. Nevertheless the estimates of the unknown parameters A, ho, hi and a% looks quite 
reasonable. I don't know yet any explanation for such success. Probably it's accidental. 

The second real solution to ( |9.4| ) is numerically quite close to the first one 

A = 0.305435 . . . 

ho = 2.44244 .. . (9.11) 
h x = -12.5161... 

The estimate ( |9.7| ) gives now z^ = —0.64821 .... It remains a question whether this solution 
is an artefact of the approximate approach or we must take it seriously and look for a proper 
interpretation. Let me note that with the second value of A ( |9.11D the "jump" between 
the asymptotic estimate (|8.1|) of h n and hi from Q9.ll ) is rather big. This makes the whole 
reasoning above very questionable. 
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10. Remarks 



The most interesting conclusion we're led by both TCS experiment and analytic con- 
siderations (in the framework of SLYM) is that the spherical partition function of a 
perturbed CFT may have quite specific analytic properties in the coupling constant 
( |3.5|) . Are these properties are typical for any perturbed CFT on a sphere? Exactly 
solvable example of free massive fermion (see Appendix A.l, where the partition func- 
tions of free massive fields are presented explicitly) shows very similar analytic struc- 
ture. Combining these two examples one might conclude that the answer is probably 
yes. However, the free massive boson (Appendix A. 2) gives an immediate counterex- 
ample. Analytic structure is more complicated, the partition function developing an 
infinite sequence of brunch points etc. Probably the simple analytic picture is charac- 
teristic for rational perturbed CFT's. It remains to be answered if integrability of the 
corresponding model in flat plays any role in this analytic structure. 

Consider the derivative z'(h) = dz(h)/dh, which is again supposed to be an entire 
function of h with the asymptotic at \z\ — > oo, —tt < axgz < 71 



z'(h) ~ exp (^Ah 5/e + ^ log/i + const + ...J 



(10.1) 



Denote its zeros h' n . Perturbative expansion ( 7.3|) gives exactly the first zero h' = 0. 
From ( |10.1| ) the large n asymptotic of h' n is read off 



1 



-h' n =(jin + ^)+0(n- 1 )) (10.2) 



6/5 



In particular, the convergent product 

~28'"- L - L r h: 



71=1 



relates z'(h) to the positions of its zeros. The last known term in the series ( |7.3| ) results 
in the sum rule for h' n 

^ 1 84 , , 

% = (10.4) 

n=l " 

Take the leading asymptotic ( |10.2|) as exact for n > 1. The sum rule becomes 

- V — — = —z 3 = 0.3862303 (10.5) 

and gives for A 

A = 0.218767... (10.6) 
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This accurate result is again remarkable, especially in view of the amazingly simple 
way it is obtained. Also, unlike more complicated system in sect. 9, this way gives a 
unique solution. 

Note, that the derivative z'(h) is related to the (unnormalized) one-point function of 
the perturbing operator 

z(h)(2R)- 2 / 5 ((f)) = -2z'(h) (10.7) 

In the TCS approach 

z(h)(2R)- 2 / 5 (0) = lim vl/tcW exp(-2f/5) (10.8) 

(no subrtactions are required for this one-point function). Apparently, for any finite- 
dimensional TCS problem this is again an entire function of h. It seems natural to 
generalize the observed analytic structure for arbitrary unnormalized one-point func- 
tion 

z(h)(2R) 2Aa ($ a ) = lim (^ ac (t) exp(2A a t) - subtractions) (10.9) 

t— »oo 

and then use it to get more information on the grounds of perturbation theory and 
analytic properties. These attempts will be reported in a separate publication |L0| . 

It would be interesting to compare both the experimental (|7.6|) and "theoretical" fl9.10|) 
estimates of the first supleading IR correction to the partition function with the pre- 
diction (|4.15|) evaluated with the use of exact form-factors of 9 [|J. 

TCS approach as described above proves to provide reasonable experimental data for 
the perturbed CFT on the sphere. It is well known however, that this method works 
well mainly for perturbed rational CFT's, and moreover, when the conformal per- 
turbative integrals are convergent. In particular, many interesting models where the 
dimension of perturbing operator is close to 1, or with asymptotically free marginal per- 
turbation are completely unaccessible by this numerical scheme. A suitable approach 
which would allow to overcome this restriction still remains to be developed. 

Nevertheless, there are still many perturbed CFT's (like the sin-Gordon model with 
sufficiently small (3) where TCS must perform reasonably and the hypotheses proposed 
in this article can be checked against the experimental data. Work in this direction is 
in progress. 
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Appendix 



A. Free fields on a sphere 
A.l. Majorana fermion 

Free massive Majorana fermion in curved geometry with conformally flat metric ( |2.4j ) is 
defined by the action 



•^•ferm 

7T 



- / d 2 x [i)di) + i)di) + ime^ 2 ^] (A.l) 



- / d 2 x [ipdip + ipdip] + — I e{x)e^ x) d 2 



x 



where m is the mass of the fermion, which is just the coupling constant in the framework of 
perturbed CFT, while the energy density 

e = 2ie^ l2 i)i\) (A.2) 

plays the role of perturbing operator. This action leads to the Dirac equations of motion 

2dip = ime^i) (A.3) 
2<9 , = —ime (p ' 2 i/) 

We're interested in the partition function Zm(R) of this field theory on the sphere of 
radius R. Obviously its derivative in m is determined by the expectation value of e. By 
symmetry (e(x)) is independent on x and 

rfl ° g 5 (it>) = -2R 2 (e) (A.4) 
am 

To evaluate (e) consider the two-point functions of the fermionic fields. For the symmetry 
reasons 

(iP{z,z)iJ;(0)) = -f(zz) (A.5) 

(■ip(z,z)ip(0)) = -ig(zz) 

Equations of motion ( |A.3| ) result in the following simple system 

(l-t)f(t) = rg(t) (A.6) 
tg'it) = rf(t) 

where r = mR and 

^ - 2 f^) (A.7) 



[1 + zz) \2RJ 
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is related to the geodesic distance s between and z. Suitable solution, which is regular at 
t = 1 (i.e., at the south pole of the sphere z = oo) and properly normalized at t = reads 

f(zz) = r(1+ir) 2 r(1 ~ ir) F (ir, -ir, 1, (1 + zz)' 1 ) (A.8) 

9(zz) = - rr( \ + J r | I ]_ 1 ) " 2r) F (1 + ir, 1 - ir, 2(1 + zz)^) 

Taking the limit zz — > one can read-off the expectation value (e(0)) at the north pole z = 

(e) = 2g(e 2 e~^)e~^/ 2 (A.9) 

777* /* 6^ 

= j (log ^ + V(l + ir) + - ir) - 2^(1) 

where ^(x) = r'(x)/r(x) and e is the invariant UV cutoff. Integrating now ( [A .4 ) and taking 
into account the conformal anomaly (|2.9|) one finds 



Z '"< fl > -espf-WlcgJL + C+Q^ ^gg^L _ (A . 10) 



z (') B i/ 6 ''v v b 2s r 2 (i + ir|i,i)r 2 (i-ir|i,i) 

where C is the Euler's constant and r 2 (x|l,l) is a particular case of the Barnes double 
gamma function [15| 

, oo 

— (A.ll) 



logr 2 (x\uJi,uj2) = -j- ^ (x + muji + nuo 2 ) 

m,n=0 



Ferimionic partition function Z m (R) is an entire function of r. It has zeros at r = ±z, 
±2z, ±3i ... of multiplicities 1, 2, 3 . . . respectively, which are of course related to the zero 
modes of the spherical Dirac operator at these values of mR. 

At small R this partition function develops as follows 

, Z ™( R ) 1, o 2, e ^ (-) fc C(2A;+ 1) 2fc+2 _ 

log ^ = -log R-r 2 log h> , V fc+2 A.12) 

& 6 & & 2# ^ fc + 1 v ; 

^0 fc=l 

The large R behavior is determined by the asymptotic expansion 



Zm (R) 



(f) 



(f) in 

where 



e X p(-4rf f « +I ^ + ss ^j...) (A.13) 



4tt V 2 2 / 
is the bulk vacuum energy of the free fermion in flat space-time and 

z® = exp(2C'(-l)) = 0.718318 . . . (A.15) 
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A. 2. Massive boson 

The action of free massive boson <ft in curved geometry reads 

Aoson = ^j d2x [(<^) 2 + ™V0 2 ] (A.16) 

and lead to the free field equation of motion 

Addcj) = m 2 e v (j) (A. 17) 

Again, the derivative in m 2 of the bosonic partition function Zm' (R) is expressed in terms 
of (0 2 ) on the sphere 

d\ogZ${R) 1 



dm 2 Arc 



[(f) eftfx = -&(<!?) (A.18) 



From (|A.17 ) it follows that the two-point function G(zz) = (<f)(z, z)(p(0)) satisfies hypergeo- 
metric differential equation 

t(l - t)G tt + (1 - 2t)G t - m 2 R 2 G = (A.19) 

where t is the same as in ( |A.7| ). Relevant solution, regular at the "south pole" z = oo and 
properly normalized at z — > 0, has the form 

s r(l/2 + is)r(l/2-is) , . , , , n . . 

G(zz) = -F (1/2 + is, 1/2 - is, 1, (1 + zz)~ l ) (A.20) 

where 

s = (m 2 i? 2 - 1/4) 1/2 (A.21) 
Considering the subleading term in z — > asymptotic of this function we find 

<0 2 ) = G (e 2 e-^°)) (A.22) 

= ~ (log ^ + ^(1/2 + is) + V(l/2 - is) - 2^(1) 

where again the UV cutoff e appears explicitly 

Integration of (|A.18|) is not that straightforward as in the fermionic case. Expectation 
value (|A.22 ) is singular at s = ±i/2, i.e, at m — 0. This is of course a manifestation of 
the zero mode of ( |A.17|) at m — 0. It leads to an extra multiplier (mi?) -1 (in addition to 
the conformal anomaly ( |2.9| )) in the small R behavior of the partition function. Normaliza- 
tion of Zrn(R) becomes ambiguous even if is fixed in some way. Below we adopt the 
normalization where Z%, (R)/Zq ~ i? 1 / 3 j (mR) at R — > 0. With this convention 

Z ™ {R) (A.23) 



4 b) # 1/3 



exp(-2C'(-l))r 2 (3/2 + is\l, l)r 2 (3/2 - zs|l, 1) 



exo (m 2 R 2 (lot ±- + C + 1^ GXP[ ^ [ ~ i))L 2[6/Z + M|1> ij 
exp (m it Jtog ^ + C + 1 j j (2 coshM)1/2 
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Bosonic partition function has brunch points of order —1/2 at s = ±z'/2, of order —3/2 
at s — ±3z/2, of order —5/2 at s = ±5z/2 etc. These singularities are caused by the 
zero modes of ( |A.17| ). In terms of the variable m 2 R 2 they turn to similar brunch points at 
m 2 R 2 = - n (n + 1), n = 0, 1, 2, . . . 
Small .R behavior reads 



~^R exp {-^r~ l log i^ + 1 )+i mR — o — m R +■■■) (a- 24 ) 



while at R — ► oo 

'(b) 



exp f-47ri? 2 ^) + 1 + 2 _\ ( A .25) 



Z m (R) _ 1/3 JV) Z' /1 _„2c(b) , ! 

7 (b) ~ m 



^ — ta^ + C-i] (A.26) 



where 



4tt V 2 2 
is the bulk vacuum energy of the massive free boson in flat, and 

zt ] = exp (1/4 - 2C'(-1)) = 1.78754. . . (A.27) 
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